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Abstract 

We illustrate some physical application of a lattice formulation of the two-dimensional 
J\f = (2, 2) supersymmetric SU(2) Yang-Mills theory with a (small) supersymmetry 
breaking scalar mass. Two aspects, power-like behavior of certain correlation func- 
tions (which implies the absence of the mass gap) and the static potential V(R) 
between probe charges in the fundamental representation, are considered. For the 
latter, for R < 1/g, we observe a linear confining potential with a finite string ten- 
sion. This confining behavior appears distinct from a theoretical conjecture that a 
probe charge in the fundamental representation is screened in two-dimensional gauge 
theory with an adjoint massless fermion, although the static potential for i? > 1/g 
has to be systematically explored to conclude real asymptotic behavior in large 
distance. 
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1 Introduction 



Recently, through the observation of a "partially conserved supercurrent rela- 
tion", we obtained [1] an affirmative numerical evidence that a lattice for- 
mulation in Ref. [2] provides a supersymmetric regularization of the two- 
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dimensional Af = (2, 2) supersymmetric Yang-Mills theory (SYM) 1 2 

S = ^Jd 2 x tr ^-F MN F MN + ^ T CT M D M <S> + # 2 } , (1) 

when one supplements to S a supersymmetry breaking scalar mass term 

Smass = ^Jd 2 x fi 2 tr {A 2 A 2 + A 3 A 3 } . (2) 

The scalar mass term was added to suppress a possible large amplitude of 
scalar fields along flat directions that may amplify 0(a) lattice artifacts to 0(1) [1]. 
In the present Letter, we illustrate some physical application of this lattice for- 
mulation for the system S + S mass . 



2 Correlation functions with power-like behavior 

Assuming the 't Hooft anomaly matching condition, in Ref . [20] , it was pointed 
out that the two-dimensional M = (2, 2) SYM has no mass gap. This as- 
pect has been numerically investigated from almost a decade ago [21,22] by 
utilizing the supersymmetric discretized light-cone formulation [23]. In this 
super-renormalizable system, it is in fact possible to determine (to all orders 
of perturbation theory) an explicit form of a correlation function between 
Noether currents, by employing anomalous Ward-Takahashi (WT) identities 
(i.e. the Kac- Moody algebra) [24]; this explicit form directly proves the above 
assertion. Here, rather than supersymmetry, continuous global (bosonic) sym- 
metries are important and the proof [24] applies even with supersymmetry 
breaking scalar mass term (2). 

The total action S + S mSuSS is invariant under the (two-dimensional) U(l)v 
transformation, \I/ — > exp{iaT 5 }\l/, and an associated Noether current (U(l)v 
current) is given by 

j, = - 2 tr {^ T CT^} . (3) 



1 For other lattice formulations of this system, see Refs. [3-9]. For recent devel- 
opments in this field of research, see Ref. [10] for a review and references cited 
in Ref. [1]. As further recent study, see Refs. [11-19]. 

2 This system can be obtained by dimensionally reducing the four-dimensional M = 
1 SYM from four to two dimensions and hence a four-dimensional notation is useful; 
Roman indices M and N run over 0, 1, 2 and 3, while Greek indices \x and v below 
run over only and 1. With the dimensional reduction, it is understood that d 2 = 
and $3 = 0. is a four-component spinor. We follow the notational convention 
in Ref. [1]. Note that the gauge coupling g has the mass dimension 1. 
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Similarly, associated with the U(1)a symmetry, \1> — > exp{aT 2 T 3 }^/ , A 2 — > 
cos{2a}A 2 — sm{2a}A 3 and A 3 — > sm{2a}A 2 + cos{2a}A 3 , there is a Noether 
current (U(1)a current), 

= j 2 tr {-i* T CT M r 2 r 3 * + 42(A 3 F M2 - A 2 F^)} . (4) 
It is then possible to show that [24], for the two-dimensional euclidean space R 2 , 

= T~(^c - !) | -7-4\9 ~~ ^uaXpXa) + c6f Jlu 5 2 (x) \ , (5) 

47T [ 7T (rj J 

to all orders of perturbation theory, where N c is the number of colors and the 
constant c is a regularization ambiguity in a divergent one-loop diagram. Thus 
the correlation function between the U(l)v current and the U(1)a current 
possesses a massless pole and this is precisely what the 't Hooft anomaly 
matching condition claims for this two-dimensional system. 

We want to confirm the power-like behavior of correlation function in Eq. (5) 
by using a lattice Monte Carlo simulation. For this, we prepared sets of uncor- 
rected configurations listed in Table 1. For simulation details, see Refs. [25,26,1] 
In the table, a denotes the lattice spacing and j3 and L are temporal and 
spatial physical sizes of our lattice, respectively. The scalar mass squared is 
[i 2 /g 2 = 0.25 for all cases. The temporal boundary condition for fermionic 
variables is antiperiodic as in Ref. [1]. For current operator (4), we discretized 
the covariant derivatives F^ 2 = D fl A 2 and = D^A^ by using the for- 
ward covariant lattice difference. Eq. (5) suggests that we should not take an 
average of the correlation function over the spatial coordinate x\ (i.e., pro- 
jection to the zero spatial momentum) because after the average, correlation 
function (5) becomes proportional to S(xq) that cannot be distinguished from 
the regularization ambiguity; we should measure the correlation function as it 
stands without the zero spatial momentum projection. 

In Fig. 1, we plotted —i(j^(x)e up j5 P (0))/2 with /i — v — along the line x\ = 0. 
We plotted also theoretical prediction (5) for M 2 with N c = 2, (3/47r 2 )l/(x ) 2 , 
by the broken line. We clearly see the power-like fall of the correlation function 
for x g < 0.7 3 instead of exponential one, although the overall amplitude is 

3 In Fig. 1, we plotted the correlation function as a function of xq, along the 
line X\ = 0. As xq moves away from the origin xq = 0, the point x approaches 
a periodic image of the origin at xq = (3 and for xo > (3/2 we expect the correla- 
tion function is power-like in the variable (3 — xq. In other words, the fact that our 
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Table 1 

Sets of uncorrelated configurations used for Figs. 1 and 2. The scalar mass squared 
is ji 2 jg 2 = 0.25 for all cases. 



lattice size 


ag 


(5g x Lg 


number of configurations 


set label 


16 x 8 


0.1768 


2.828 x 1.414 


400 


I 


20 x 10 


0.1414 


2.828 x 1.414 


800 


II 


24 x 12 


0.1179 


2.828 x 1.414 


400 


III 


20 x 16 


0.1414 


2.828 x 2.263 


400 


IV 




Fig. 1. The correlation function —i{j ll {x)e vp j^p{ l d))/{2g 2 ) with fi = u = along the 
line x\ = 0, for the configuration sets in Table 1. The broken line is theoretical 
prediction (5) for M 2 . 

somewhat larger than the theoretical expectation for R 2 . From the behavior in 
the figure, we think that this discrepancy in the overall amplitude is caused by 
a finite lattice spacing and volume. In particular, comparison between set II 
(indicated by x) and set IV (indicated by O) shows that the finite size effect 
is rather large (note that these two sets differ only in the spatial physical 



finite-size lattice is topologically T 2 but not R 2 cannot be neglected for xo > (3/2. 
We thus do not expect the power-like fall (that is expected for 1R 2 ) for x$g > 1 and 
actually the plot blows up for xog > 1 (in our simulation, fig = 2.828). This remark 
is applied also to Fig. 2, in which the antiperiodic boundary condition for fermionic 
fields implies "blow-down" for x^g > 1. 
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size L). We thus expect that the theoretical prediction for R 2 is eventually 
reproduced in the limit, a — > and /3, L — > oo, although we do not carry out 
a systematic study on this limit. 

What is the implication of the above observation? It indicates that our target 
theory, the two-dimensional M = (2, 2) 577(2) SYM with a scalar mass term, 
is realized in the continuum limit of the present lattice model. In particular, 
in deriving Eq. (5), one assumes that the U (l)y and U (1)a currents j M and 
individually conserve [24]. 4 One assumes U(l)v and U(1)a symmetries in this 
sense. In the present lattice formulation [2], the U(l)y symmetry is explicitly 
broken for finite lattice spacings. The above observation hence indicates that 
the U(l)y symmetry is fairly restored with present lattice spacings. (This 
symmetry will eventually be restored in the continuum limit [1].) 

Now, if the system were supersymmetric, and if supersymmetry is not sponta- 
neously broken, there would exist a massless fermionic state corresponding to 
the massless bosonic state appearing in Eq. (5) as an intermediate state. We 
expect that this fermionic state produces a massless pole in the correlation 
functions 

((fy)»(aO(/i/)i(0)) (i = 1, 2, 3, 4; no sum over i), (6) 
where i denotes the spinor index and 

s M = -\cr M r N r p tr {f mn ^} , (7) 
fn = ^ (r 2 tr{A 2 v} + r 3 tr{ }) . (8) 

In the above, s M is the supercurrent associated with the supersymmetry of S, 
5 A M = ie T CT M y, 5^ = ^F MN Y M T N e + iHT 5 e, and SH = -ie T CT 5 T M D M ^, 
and is a lowest-dimensional fermionic spinor- vector (considered in Ref. [1]). 
Eq. (6) with % — 1, 2, 3, and 4 are precisely four correlation functions studied 
in Eq. (11) of Ref. [1] and, as noted there, these four functions are identical to 
each other in the continuum theory. Our expectation that Eq. (6) possesses a 
massless pole stems from a supersymmetric WT identity for a vanishing scalar 
mass squared, /i 2 = 

^ i=i 

= ~\ MxfaMO)) - L(x)e±tr{A 3 (0)F p2 (0) - A 2 (0)F p3 (0)}\ , (9) 



4 This assumption fails for example in the massless Schwinger model, in which the 
U(1)a current suffers from the axial anomaly; note that the massless Schwinger 
model has a mass gap. 
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which follows from 5(j^{x) /J(0)) = 0, where 5 is the global super transforma- 
tion; this relation holds under the assumptions that the boundary condition 
is consistent with supersymmetry and supersymmetry is not spontaneously 
broken. (In deriving Eq. (9), we have used also the equation of motion of the 
auxiliary field, H = H — iF 01 = 0). In the right-hand side of Eq. (9), the mass- 
less pole in the first term (recall Eq. (5)) cannot be cancelled by the second 
term, because the latter is 0(g 2 ) as one can easily see. 5 

Even if the supersymmetry breaking owing to lattice regularization disappears 
in the continuum limit [1], our present system is not supersymmetric because 
there is scalar mass term (2) and we used the antiperiodic temporal bound- 
ary condition for fermions. These will give additional contribution to Eq. (9). 
In Fig. 2, we plotted correlation functions (6) along the line x± — for set IV 
in Table I. (For the parameters of this configuration set, naively-expected or- 
der of magnitude of supersymmetry breaking caused by above factors would 
be ~ /i = 0.5(7 an d ~ 1//3 — 0.3536(7, respectively.) For 0.2 < x$g < 1.0, 
for all i, the power-like behavior expected from supersymmetric WT iden- 
tity (9) combined with Eq. (5) is fairly observed. Somewhat surprisingly, we 
do not see a significant effect of the supersymmetry breaking and it appears 
that the fermionic intermediate state is approximately massless as expected 
from approximate supersymmetry. This result is consistent with the conclusion 
of Ref. [1] that the supersymmetry breaking owing to the lattice regularization 
disappears in the continuum limit. 



3 Potential energy between probe charges in the fundamental rep- 
resentation 

Contrary to naive intuition, it is believed that a probe charge in the fun- 
damental representation is screened by dynamical adjoint massless fermions 
in the two-dimensional SU(N C ) QCD [27,28]. This phenomenon is analo- 
gous to the screening of a fractional charge in the massless Schwinger model 
with an integer-charged fermion and is believed to occur also in the two- 
dimensional Af = (1,1) SYM, despite the presence of a scalar field and a 
Yukawa interaction in the latter [27,29] (see also Ref. [30]). As a generaliza- 
tion of these, in Refs. [29,30], it was claimed that this screening persists in 
any two-dimensional (supersymmetric and non-supersymmetric) gauge the- 

5 In fact, a one-loop calculation of Eq. (6) on the basis of Eq. (1) coincides with 
expression (5), possibly with a different regularization-dependent constant c. Note 
also that even if correlation functions (6) possess a massless pole of structure (5), this 
does not necessarily imply the existence of the Nambu-Goldstone fermion associated 
with the spontaneous supersymmetry breaking. This is because structure (5) shows 
that the Fourier transform of d lx ((s lx )i(x)(f 1/ )i(0)) vanishes at zero momentum. 
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Fig. 2. The correlation function ((s fJi )i(x)(f u )i(0))/g 2 with \i = v = along the 
line x\ = for set IV in Table 1. The broken line is (3/47r 2 )l/(xo) 2 , the same 
function plotted in Fig. 1. 

ory with adjoint massless fermions, although an explicit proof was not given 
there. In our present system, the masslessness of the gaugino is ensured by the 
global U(l) A and U(l) v symmetries and, hence, it is of interest to study the 
static potential energy between probe charges in the fundamental representa- 
tion. If the expected screening occurs, the static potential would approach a 
constant for large distance (i.e., the Wilson loop obeys the perimeter law). 

We thus measure the expectation value of the Wilson loop, 

W(T,R) = (^tT^J[U^y (10) 

where C denotes a rectangular loop of a physical size T x R and link vari- 
ables Ue belong to the fundamental representation of the gauge group SU(2). 
For this average, we prepared uncorrelated configurations listed in Table 2 
(the scalar mass squared is fi 2 /g 2 = 0.25 for all cases). Theoretically, the 
static potential V(R) is defined by the asymptotic form in T — > oo: 

-\n{W(T,R)} = V(R)T + c(R). (11) 

Practically, with finite-size lattices, we made a linear x 2_n t of — ln{VF(T, R)} 
with respect to T in a finite range T min < T < (3/2 for each R and regarded 
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Table 2 

Sets of uncorrelated configurations used for Figs. 3 and 4. The scalar mass squared 
is ji 2 jg 2 = 0.25 for all cases. 



lattice size 


ag 


09 x Lg 


number of configurations 


set label 


20 x 10 


0.2 


4x2 


800 


V 


20 x 10 


0.1414 


2.828 x 1.414 


800 


VI 


20 x 16 


0.1414 


2.828 x 2.263 


400 


VII 


30 x 10 


0.1414 


4.243 x 1.414 


800 


VIII 



the slope as V(R); obviously (5/2 is a maximal temporal size of the Wilson 
loop that is physically meaningful. We determined the lower end of the fit T min 
such that the fitting range becomes as wide as possible insofar as % 2 /dof of 
the fit does not exceed unity. We had T min = a - 5a. A typical result of this 
linear fit is depicted in Fig. 3 for the case of set VI in Table 2. 




k' : - ' 1 1 1 1 1 1— 

0.2 0.4 0.6 0.8 1 1.2 1.4 

Tg 



Fig. 3. The linear x 2 -ht of - ln{W(T, R)} for set VI in Table 2. 

In Fig. 4, we plotted V(R) for R < L/2 determined in this way for various 
lattice spacings and lattice sizes (Table 2). The error in the figure was de- 
termined by the range of a slope of the linear fit that corresponds to a unit 
variation of x 2 - 

Now in Fig. 4, all points are almost on a common line, although lattice spac- 
ings are different (ag = 0.2 for + and ag = 0.1414 for x, □ and Q). This fact 
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0.3 



0.25 - 



set V i — i — i 
VI 

VII r-B—: 

VIII : O--! 



m 



0.2 - 



0.15 - 




# 



0.1 - 



+ 



0.05 - 



H 











0.2 



0.4 



0.6 

Rg 



0.8 



1 



1.2 



Fig. 4. V(R)/g determined by the linear x 2 -fh described in the text. See Table 2 
for the label of configuration sets, /j, 2 j g 2 = 0.25. 

indicates that the result in Fig. 4 can roughly be regarded as that in the con- 
tinuum limit. Similarly, since physical lattice sizes of each configuration set are 
rather different (for example, Lg = 1.414 for x and Lg = 2.263 for □), there 
appears almost no significant finite-size effect. 6 Therefore, at least for Rg < 1, 
we could conclude that the static potential V(R) is linear (i.e., the Coulomb 
potential in two dimensions) with a finite string tension o ~ 0.25g 2 . This ap- 
pears to be distinct from a theoretical conjecture in Refs. [29,30]. 7 However, 
to conclude whether a probe charge is really confined or screened, the static 
potential for Rg > 1 has to be systematically explored; we reserve this as a 
future project. 

It is also of interest to see how the behavior in Fig. 4 changes as a function 
of the scalar mass. For example, in the limit fi 2 /g 2 — > oo, the scalar fields 
will completely decouple 8 and our system would become the two-dimensional 



6 The discrepancy between □ and x at Rg = 0.5657 could be explained by the fact 
that this value of Rg is comparable with the spatial lattice size for set VI. Note that 
for smaller Rg they have less discrepancies. 

7 A possible confutation is that the gaugino is not strictly massless in our simulation 
because of the antiperiodic temporal boundary condition. This point seems irrele- 
vant, however, because the behavior in Fig. 4 appears insensitive to the temporal 
size (3 of our lattice. Compare, for example, set VI and set VIII. 

8 A unique UV divergent diagram that contains scalar loops is a one-loop scalar self- 
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Table 3 

Sets of uncorrelated configurations used for Fig. 5. 





lattice size 


ag 


09 x Lg 


number of configurations 


1.69 


20 x 10 


0.1414 


2.828 x 1.1414 


800 


0.04 


20 x 10 


0.1414 


2.828 x 1.1414 


800 



577(2) QCD with an adjoint massless fermion. On the other hand, the limit 
A* 2 /*? 2 ~~ * would provide a possible definition of the two-dimensional M = 
(2, 2) SYM. It is believed that the screening occurs in both theories, as al- 
ready noted. To have a rough idea on this issue, we carried out a preliminary 
experiment by using sets of configurations listed in Table 3. The results are 
summarized in Fig. 5. For both [i 2 /g 2 = 1.69 and fi 2 /g 2 = 0.04, we still see 



0.35 



0.3 - 



0.25 - 



0.2 h 
0.15 

0.1 - 
0.05 - 




0.1 



□ 



+ 
x 



* 

X 



|x 2 /g 2 =1.69 i — i- 

0.25 h- x- 

0.04 :— Q- 
i i 



0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



Rg 



Fig. 5. V{R)/g determined by linear x 2 -fit described in the text. We used the con- 
figuration sets in Table 3 and, for fi 2 /g 2 = 0.25, set VI of Table 2. Eq. (12) with 
N c = 2 is plotted by the broken line. 



a linear potential, although the string tension appears somewhat smaller for 



energy. This contributes to simply shift the tree-level mass /i 2 by ~ g 2 ln{fi 2 /A 2 }, 
where A is the UV cutoff, and does not affect a complete decoupling of the scalar 
fields in the limit jJ? /g 2 — ► oo. 
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smaller /i 2 /g 2 . 9 In the figure, just for reference, we also plotted the function 





V(R) = K -±g 1 - exp -J^gR (12 



with N c = 2, that is given by a semi-classical analysis of a bosonized version 
of the two-dimensional massless QCD [27]. Strictly speaking, the overall pro- 
portionality constant is not determined by this analysis and we have chosen 
it as above without any special reason. 

Fig. 5 is simply a result with a single lattice spacing and a single lattice size. 
It is thus not clear what is the real behavior in the continuum and the large 
volume limits. Our result is still preliminary and a further detailed numerical 
study is needed. 



4 Conclusion 



In this Letter, we illustrated some numerical use of the lattice formulation [2] of 
the two-dimensional M = (2,2) SYM with a (small) supersymmetry breaking 
scalar mass. Two physical problems were considered. For the first one (Sec. 2), 
our Monte Carlo result fairly reproduced theoretical prediction on the basis of 
global symmetries and (approximate) supersymmetry in the continuum theory. 
For the second one (Sec. 3), our result for the static potential V(R) did not 
exhibit the screening behavior that theoretically anticipated. However, since 
our result of V(R) was limited for Rg < 1, it is desirable to carry out a further 
systematic study by using finer and larger lattices. 

We would like to thank Koji Hashimoto and Daisuke Kadoh for helpful dis- 
cussions. We thank also the authors of the FermiQCD/MDP [31,32] and of a 
Remez algorithm code [33] for making their codes available. Our numerical re- 
sults were obtained using the RIKEN Super Combined Cluster (RSCC). I. K. 
is supported by the Special Postdoctoral Researchers Program at RIKEN. The 
work of H. S. is supported in part by a Grant-in-Aid for Scientific Research, 
18540305. 



9 Consider the case pL 2 /g 2 = 0.04 in Fig 5. For this case, the Compton wavelength 
of a (free) scalar particle is l//x = 5.0/g and this is several times longer than the 
physical lattice size. Thus in this case the scalar field could effectively be regarded 
as massless and the points □ might be regarded as those for the two-dimensional 
M = (2, 2) SYM. 
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